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1 Introduction 

co . 

Systems of elements with motivated behavior (systems with motivation), e.g., 
fish and bird swarms, car ensembles on highways, stock markets, etc. often dis- 
play noise- induced phase transitions (for a review see Ref. [3]). The ability of 
noise to induce phase transitions is now well established (see, e.g., Refs [1,2]). 
■ However, in systems with motivation there is a special mechanism endowing 

£5 ' the corresponding noise-induced phase transitions with distinctive properties. 

For example, people as elements of a certain system cannot individually 
control all the governing parameters. Therefore one chooses a few crucial pa- 
rameters and focuses on them the main attention. When the equilibrium with 
respect to these crucial parameters is attained the human activity slows down 
retarding, in turn, the system dynamics as a whole. For example, in driv- 
ing a car the control over the relative velocity v is of prime importance in 
comparison with the correction of the headway distance x. So, under normal 
conditions a driver should eliminate the relative velocity between her car and 
a car ahead first and only then correct the headway. 
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Fig. 1. Evolution of the distribution function V(x,v) (shown by level contours) as 
the parameter A decreases. The parameters used are a = 0.1 and e = 0.1. The lower 
right window depicts only one maximum of the distribution function. 



These speculations lead us to the concept of dynamical traps, a certain 
"low" dimensional region in the phase space where the main kinetic coeffi- 
cients specifying the characteristic time scales of the system dynamics become 
sufficiently small in comparison with their values outside the trap region [4, 5]. 
The present paper analyzes the effect of noise on such a system and demon- 
strates that additive noise in a system with dynamical traps is able to give 
rise to new phases. 



2 Noised oscillatory system with dynamical traps 

By way of example, the following dimensionless system typically used to de- 
scribe the oscillatory dynamics is considered: 

^=v, ^ = -Q{v) [x + av] + e£(t) , (1) 

where a is the damping decrement and the term e£(i) is a random Langevin 
"force" of intensity e proportional to the white noise with unit amplitude. 
The function f2(v) describes the dynamical trap effect arising in the vicinity 
of v = 0. For this function, the following simple Ansatz 

v 2 + A 2 

v z + 1 

is used. In the chosen scales the thickness of the trap region is equal to unity 
and the parameter A < 1 measures the trapping efficacy. When A = 1 the 
dynamical trap effect is ignorable, for A = it is most effective. 

System (1), (2) was analyzed numerically using the algorithm described 
in [6]. Figure 1 shows the distribution function V(x,v) of the system on the 
phase plane {x, v}, depending on the parameter A. As can be seen this system 
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Fig. 2. A typical fragment of the system path going through the trap region. The 
parameters a = 0.1, e = 0.1, and A = 0.01 were used in numerical simulations in 
order to make the trap effect more pronounced. 

undergoes a second order phase transition manifesting itself in the change of 
the shape of the phase space density V(x, v) from unimodal to bimodal as the 
trap parameter A decreases. 

Mechanism of the phase transition 

Figure 2 illustrates the mechanism of this phase transition. It depicts a typical 
fragment of the system motion through the trap region for A < 1. When the 
path goes into the trap region Q t (\v\ <C 1) the regular "force" Q(v)(x + av) 
is depressed. So inside this region the system dynamics is mainly random due 
to the remaining weak Langevin "force" e£(£). Crucial is the fact that the 
boundaries d+Qt (where v ~ 1) and d-Qt (where v ~ — 1) are not identical 
in properties with respect to the system motion. At the boundary <9+Qt the 
regular "force" leads the system inwards the trap region Q t , whereas at the 
boundary d-Q t it causes the system to leave the region Q t . Outside the trap 
region Qt the regular "force" is dominant. Thereby, from the standpoint of 
the system motion inside the region Qt, the boundary d+Qt is "reflecting" 
whereas the boundary cL Qt is "absorbing" . 

As a result the distribution of the residence time at different points of 
the region Q t should be asymmetric as shown in Fig. 1 (it is most clear in 
the lower right window). Therefore, during location inside the trap region the 
mean velocity must be positive, causing the system to go away from the origin. 
Outside the trap region the system motion is damping. So, when the former 
effect becomes sufficiently strong the distribution function V(r\, u) becomes 
bimodal. 
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n s (Vi-i,Vi,V i+1 ) = — - , (5) 



3 Chain of oscillators 

A similar noise-induced phase transition for an ensemble of oscillators with 
dynamical traps is analyzed. The following one-dimensional model is consid- 
ered, TV balls can move along x-axis interacting with the nearest neighbors, 
which is described by the equations for i = 1, 2, ... N — 1 

dx ' dv ' 

-jf =v i> -j£ = -fts(vi-i,Vi,Vi +1 )[r]i+<T'&i]+e£i(t). (3) 

Here X{ is the coordinate of ball i, v% is its velocity, the variables m and i?j are 
given by the expression: 

rn =Xi - \{xi-i + x i+1 ), §i = Vi - 5(^-1 + v i+ i) , (4) 

and {£»(t)} is the collection of white noise sources of unit amplitude and being 
mutually independent. The damping decrement a and the noise intensity e are 
assumed to be the same for all the oscillators. The boundary balls (i = and 
i = N) are let to be fixed to prevent the ball system from moving as a whole. 
The function f2 s (vi-i,Vi,Vi+i) measuring the trap effect due to the nearest 
neighbor interaction is given by the Ansatz 

(v^! - V t ) 2 + (Vt - v l+1 ) 2 + A 2 
(vi-i - Vi) 2 + (vi - v i+ i) 2 + 1 

where the parameter A has the same meaning as previously, it measures the 
intensity of trapping. The balls are assumed to be either mutually permeable 
or impermeable. In the latter case the absolutely elastic collision approxima- 
tion is used. 

The system of equations (3)-(5) was analyzed numerically. Below, the 
results are presented for the impermeable balls. Again integration of the 
stochastic differential equations was performed with the algorithms described 
in [6] . We analyzed an ensemble of 500 balls initially spaced 5 units apart and 
counted the number of balls falling into 100 identical intervals covering the 
system location region in the corresponding space. 

Figure 3 actually depicts the ball distribution along the spatial axis and 
the velocity space at a fixed moment of time. As seen, the trap effect leads 
to the formation of a sufficiently inhomogeneous spatial distribution of balls 
whereas for the same system but without trapping, A = 1, the balls arc 
actually uniformly distributed over the space region. In the velocity space the 
trap effect causes the balls to spread over a much wider domain. 

To make the trap effect in the velocity space more clear we averaged the 
velocity distribution over 500 time units. The result is presented in Fig. 4. 
As seen, the strong trap effect causes a nonanalytic behavior of the velocity 
distribution at zero value, it takes the form of a cusp. 

Finally, Fig. 5 illustrates the time pattern of the system dynamics. It is 
clear that the shown spatial structure is characterized by a long life time, 
because, we recall, the period of individual ball oscillations without traps is 
about 2ir. 
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Fig. 3. The ball distribution along the spatial axis and the velocity distribution at 
a fixed moment of time. These distributions are represented in terms of number of 
balls falling into measuring interval. The measuring intervals of number 100 cover 
the whole region of the system location in the corresponding space. The thick grey 
line and the boundary of dashed region match the cases where the trap effect is 
absent, A = 1, and strong, A = 0.1, respectively. The other parameters used are 
a = 0.1 and e = 0.1. 
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Fig. 4. The time averaged velocity distribution presented in the same way as in 
Fig. 3. Here, however, the individual length of the measuring interval was decreased 
so their total number was 250. The parameters used are A = 0.1, a = 0.1 and 
e = 0.1, the averaging time interval is 500 unit. The left and right windows depict 
the same data in linear and logarithmic scales, respectively. 



4 Conclusion 

A new type of noise-induced phase transitions in systems of elements with 
motivated behavior is considered. Dynamics of such systems, we think, should 
exhibit a number of anomalies due to dynamical traps. The dynamical traps 
form a low dimensional region in the phase space where the kinetic coefficients 
become sufficiently small and, as a result, the system spends a long time 
in it. The cause of the dynamical traps is, e.g., the inability of people or 
animals to control all the system governing parameters. So they have to focus 
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Fig. 5. A fragment of the system time pattern. Black regions correspond to increased 
number of balls located in them. 

the main attention on a few crucial ones and the intensity of their activity 
decreases when the system attains a local quasi-equilibrium with respect to 
these parameters. 

To illustrate this effect a simple oscillatory system {x,v — x} is studied 
when the trap region is located in the vicinity of the rr-axis and without noise 
the stationary point {x = 0, v = 0} is absolutely stable. For this system as 
shown numerically an additive white noise can cause the phase-space density 
to take a bimodal shape. For the chain of such oscillators the trap effect gives 
rise to a substantially nonuniform spatial distribution and leads to nonanalytic 
behavior of the velocity distribution near zero value. 

It should be underlined that a possible phase state that could be ascribed 
in this case to a maximum of the distribution function in the phase space does 
not match any stationary point of the "regular" or "random" forces. 
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